Invariant imbedding theory of mode conversion in inhomogeneous plasmas: I. Exact 
calculation of the mode conversion coefficient in cold, unmagnetized plasmas 
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This is the first of a series of papers devoted to the development of the invariant imbedding theory 
of mode conversion in inhomogeneous plasmas. A new version of the invariant imbedding theory 
of wave propagation in inhomogeneous media allows one to solve a wide variety of coupled wave 
equations exactly and efficiently, even in the cases where the material parameters change discon- 
tinuously at the boundaries and inside the inhomogeneous medium. In this paper, the invariant 
imbedding method is applied to the mode conversion of the simplest kind, that is the conversion 
of p-polarized electromagnetic waves into electrostatic modes in cold, unmagnetized plasmas. The 
mode conversion coefficient and the field distribution are calculated exactly for linear and parabolic 
plasma density profiles and compared quantitatively with previous results. 



I. INTRODUCTION 



The conversion of one type of wave mode into 
another type of mode at resonance points in inho- 
mogeneous plasmas and the associated irreversible 
transfer of wave energy to the mode conversion 
region arc important and extensively-studied phe- 
nomena in both space and laboratory plasmasiSiMi^ 
There exists a vast literature on a large number of 
mode conversion phenomena involving many different 
wave modes, which we will not attempt to review 

l ier6| 6 i 7 i 8 4 9 4 10 i ll i 12 i 13 4 14 i 15 4 l^ 

Since mode conversion is associated with the singularity 
of wave functions, theoretical studies encounter great 
difficulties and often adopt some kind of approximate 
methods such as the WKB method. In complicated 
cases where plasma waves are strongly coupled, these 
approximate methods are often unable to produce 
sufficiently accurate results. Therefore development of a 
theoretical method that allows exact solutions to mode 
conversion problems will be of great value in the study 
of inhomogeneous plasmas. 

In this paper and a series of companion papers to 
be published later, we will develop an exact theory of 
mode conversion based on a new version of the invari- 
ant imbedding method (IIM) for wave propagation in 
arbitrarily-inhomogeneous stratified media. The main 
idea of the IIM is to transform the original bound- 
ary value problem of wave equations, which are second- 
order differential equations, to an initial value problem 
of coupled first-order ordinary differential equations for 
the reflection and transmission coefficients and the field 
amplitudesi 29 ' 30 i 31 ' 32 i 33 ' 34 i 35 These equations are called 
the invariant imbedding equations and the independent 
variable, which is usually the thickness of inhomogeneous 
media, is called the imbedding parameter. This transfor- 
mation makes the numerical solution of wave equations 
much easier. 

Our IIM utilizes a rigorous integral representation of 



wave equations and differs substantially from the better- 
known IIM presented in Ref. 29 l 17 i 18 i 33 i 34 i 35 More specif- 
ically, our method gives the invariant imbedding equa- 
tions with nonsingular coefficients even in the cases where 
the material parameters change discontinuously at some 
discrete points inside the inhomogeneous medium, un- 
like the imbedding equations derived in Ref. 29. When 
the inhomogeneity is one-dimensional, which is the case 
in stratified media, our method can be used to obtain 
exact solutions for the reflection and transmission coef- 
ficients and the electric and magnetic field amplitudes 
inside arbitrarily-inhomogeneous media. When the inho- 
mogeneity is random, it can be used to obtain the exact 
disorder-averaged reflection and transmission coefficients 
and field amplitudes. 33 Our method can also be used to 
study the propagation of several coupled waves in both 
linear and nonlinear media in an exact manner^ 

In this paper, we apply our IIM to the mode conver- 
sion of the simplest kind, that is the conversion and res- 
onant absorption of obliquely-incident p-polarized elec- 
tromagnetic waves into electrostatic modes in cold, un- 
magnetized plasmasi 19 i 20 i 21 i 22 i 23 i 24 i 25 i 26 i 27 i 28 When a p 
wave of frequency to enters a stratified plasma with a 
monotonically-increasing density profile at an incident 
angle 9, it propagates first to the cutoff point where the 
local plasma frequency u p equals lo cos 9 and is partially 
reflected. Some fraction of the wave then tunnels to the 
mode conversion point where u> p is equal to uj and is con- 
verted to electrostatic modes. Since the generation of 
electrostatic modes requires the electric field component 
in the direction of inhomogeneity, only obliquely-incident 
p waves can produce mode conversion in unmagnetized 
plasmas. 

Mode conversion in unmagnetized plasmas has been 
studied extensively over many decades and a detailed re- 
view of the literature can be found in Ref. 25. In this 
paper, we restrict our attention to cold plasmas with 
zero temperature and postpone the discussion of finite 
temperature effects to future publications. Though our 
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method can be applied to stratified plasmas with arbi- 
trary density profiles very easily, we consider only the 
cases with linear and parabolic density profiles in this 
paper. One of the main reasons why we revisit this 
well-studied problem is to establish the accuracy and ef- 
ficiency of our method before applying it to more com- 
plicated mode conversion problems. 

Among a large number of references on the mode con- 
version in unmagnetized plasmas with a linear density 
profile in half- or entire space, we mention just two pa- 
pers. The readers may refer to Ref. 25 for a more com- 
plete list of references. Forslund et al. performed a nu- 
merical simulation of warm, unmagnetized plasmas at 
low temperatures; 2 - The mode conversion coefficient and 
the field distribution at the lowest temperature they have 
studied are generally considered to be an accurate solu- 
tion of the mode conversion problem in cold, unmagne- 
tized plasmas. 

Hinkel-Lipsker et al. derived analytical formulas 
for the mode conversion coefficient in the cases where 
the density profile is linear and parabolic in entire 
g p ace| 25 i 26 i 27 i 28 Their result in the linear case agrees 
pretty well with Forslund et aZ.'s result, but there is some 
numerical discrepancy^ In the parabolic case, they con- 
sidered the case where the wave frequency is fixed at the 
value of the local plasma frequency at the peak of the 
parabolic density profile separately from more general 
cases where the wave frequency is lower than the peak 
plasma frequencyi^Si2L^ 

We calculate the mode conversion coefficient and the 
magnetic field distribution for precisely the same models 
as those considered in Refs. 24—28. We find a very good 
agreement with the result of Ref. 24 in the linear case, 
though we believe that the numerical accuracy of our 
result is better. In the parabolic case, we find substantial 
differences between our results and those of Refs. 26—28. 
To the best of our knowledge, the results presented in 
this paper represent the first exact solution of the mode 
conversion problem in cold, unmagnetized plasmas with 
linear and parabolic density profiles. 

In Sec. [HI we introduce the wave equation and the 
plasma density profiles. In Sec. IIIII we present the IIM 
and the invariant imbedding equations for the reflection 
and transmission coefficients and the field amplitudes. 
In Secs. lIVl andlVl the results on the case with a linear 
density profile in entire or half- space are presented. In 
Sec. I VII the results on the case with a parabolic density 
profile in entire space are presented. We conclude the 
paper in Sec. IVIII 



tivity e varies only in one direction in space. We take 
this direction as the z axis and assume that the inhomo- 
geneous plasma lies in < z < L and the wave prop- 
agates in the xz plane. Since the plasma is uniform in 
the x direction, the x component of the wave vector, q, 
is a constant and the dependence on x is taken as being 
through a factor e tqx . When q is nonzero, the wave is 
said to pass through the plasma obliquely. 

For q 7^ 0, we need to distinguish two independent 
cases of polarization. In the first case, the electric field 
vector is perpendicular to the xz plane and the magnetic 
field vector lies in that plane. This type of wave is known 
as an s (or TE) wave. In the second case, the magnetic 
field vector is perpendicular to the xz plane and the elec- 
tric field vector lies in that plane. Then the complex 
amplitude of the magnetic field, B — B(z), satisfies 



d 2 B 
~d^ 



1 de dB 

e(z) dz dz 



+ [k 2 e{z) -q 2 ]B = 0. (1) 



This type of wave is called p (or TM) wave. In Eq. QJ, 
we have assumed that the magnetic permeability of the 
plasma is equal to 1. It is well-known that if the incident 
wave is s-polarized, mode conversion does not occur in 
unmagnetized plasmas. In this paper, we will consider 
only the p wave case. 

It is straightforward to derive an analytical expression 
for e in a cold, unmagnetized electron plasma: 



e(z) 



Aire' 



uj (lu + i"f) ' 



n(z), 



(2) 



where 7 is the collision frequency and oj p is the electron 
plasma frequency, e and m are the electron charge and 
mass, respectively. We assume that the wave is incident 
from region I where z > L and transmitted to region 
II where z < 0. The dielectric permittivity is equal to 
t\ in region I and ti in region II. The phenomenological 
collision frequency 7 is assumed to be zero in regions I 
and II. 

Our method can be applied to the general case where 
the electron density n is an arbitrary function of z. In 
this paper, we restrict our consideration to three special 
cases. In Model A, the density depends on z linearly in 
entire space. In order to treat this case, we assume that 
the density n(z) is given by 



n A {z) = n 1 



A 



(3) 



II. MODEL 

We consider a plane, monochromatic and linearly- 
polarized electromagnetic wave of frequency to and vac- 
uum wave number kg = oj/c, incident from a homo- 
geneous region on a stratified inhomogeneous plasma, 
where the electron density n and the dielectric permit- 



for < z < L and take the L — > 00 limit numerically. 
The constant A is the scale length for the linear density 
profile. We point out that for z > A + L/2, the den- 
sity takes unphysical negative values. When L — > 00, 
ha{z) varies from —00 at z = L to +00 at z = 0. 
We make a further simplification by assuming that the 

frequency is fixed at the value of the local plasma fre- 

1/2 

quency at z = L/2, which is (47rnoe 2 /m) , and 7 <C uj. 
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We make these inessential assumptions just to make our 
model look exactly the same as the corresponding model 
in Ref. 25. Then the dielectric permittivity takes the 
form 



XL 
2 A 



2 A 



IT) 



if z > L 
if < z < L 
if z < 



(4) 



where the dimensionless damping parameter rj is propor- 
tional to 7 and will be sent to zero numerically We note 
that in the absence of damping, the dielectric permittiv- 
ity vanishes at z = L/2, where the mode conversion takes 
place. 

In Model B, we assume that the wave is incident from a 
vacuum on an inhomogeneous plasma, where the electron 
density depends on z linearly in half-space. The density 
is given by 



n B {z) = n Q 



L 



A 



(5) 



for < z < L. In the L — ► oo limit, tib(z) varies from 
at z — L to +oo at z = 0. Similarly to the Model 
A case, we make an assumption that the frequency is 
fixed at the value of the local plasma frequency at z = 

1/2 

L — A, (47rnoe 2 /m) , and 7 <C ui. Then the dielectric 
permittivity has the form 



1 



e B {z) = { 1 



1 - £ 

1 A 



-L 
A 



11] 



if z > L 
if < z < L 
if z < 



(6) 



The mode conversion occurs at z = L — A in this case. 

In Model C, we assume that the density depends on z 
parabolically in entire space and is given by 



n c {z) = n 



A 



(J) 



for < z < L and the L — > 00 limit will be taken nu- 
merically, no is the electron density at the top of the 
parabolic density profile and A is the scale length. As- 
suming that 7 oj, the dielectric permittivity is given 
by 



v2 Ay 



*c(z) = \ (<5 + l)(^i)' 



it] 



'Hi 
,2 Ay 



if z > L 
if < z < L 
if z < 



(8) 

where S — Aim^e 2 / '(raw 2 ) — 1. If the frequency is fixed at 
the value of the plasma frequency at the top of the density 
profile, S is equal to zero. When the wave frequency 
is shifted downward from the peak plasma frequency, 5 
takes a positive value. Then the mode conversion occurs 
at z = L/2 ±[8/ (<5 + l)] 1/2 A. 

We point out that strictly speaking, Models A and C 
make sense only in the limit where the damping parame- 
ter rj goes to zero. For any finite value of 77, the incoming 
wave of finite amplitude will decay away before reaching 
the mode conversion point. When we present our numer- 
ical result for 77 values which are not sufficiently small, we 
will always specify the thickness of the plasma, L, used 
in the calculation. 



III. INVARIANT IMBEDDING EQUATIONS 



We consider a plane 
H(x,z) = H{z)e iqx = 



p wave of unit magnitude 



Jp(L-z)+ lq x where 



P 



le\ ko cos 9 and q = fco sin8, incident on the plasma 
from the right (that is, the region where z > L). 9 is 
called the angle of incidence. The quantities of main 
interest are the complex reflection and transmission co- 
efficients, r = r(L) and t = t(L), defined by the wave 
functions outside the medium: 



H{x,z) 



e ip(L-z)+iqx _|_ r r£\ e ip(z-L)+iqx ^ Z >L 
t(L)e- ip ' z+lqx , Z<0 
I 



(9) 



where p' (= ^/e^ kcos9') is a positive real constant de- 
fined by the dielectric permittivity in region II, 62 (> 0), 
and the angle that outgoing waves make with the nega- 
tive z-axis, 9' . If £2 < 0, the wave is not transmitted to 



region II and the transmission coefficient needs not to be 
considered. 

Using the IIM, we derive exact differential equations 
satisfied by r and ti 17 i 34 i 35 



dr(l) dl) i 

= 2iy/ei k cos9 r(l) — -V £ i k cos9 

dl €1 2 

dt(l) e(l) i 

— — = i^/elk cos9 t(l) - -y/ei k cos9 

at e\ I 



e(0 
e(l) 



61 tan 2 1 



- 1 



e(l) 
e(l) 



1 - tan 2 1 



[1+Kor, 

l + r(l)]t(l). 



(10) 



4 



These equations are supplemented with the initial condi- 
tions for r and t, which are obtained using the well-known 
Fresnel formulas: 



r(Q) 
f(Q) 



e 2 Jei cos 



0-er/ 



e 2 - (\ sin 



C2y/ei cos 8 + eiV£2 — ei sin 2 
2e 2v / ei cos0 



e 2x /£i cos 



e 2 - ei sin 



2 , 



(11) 



When (e 2 — ei sin 2 9) is negative, the square roots appear- 
ing in Eq. (fTTfl have to be replaced by i\ 



|e 2 - ex sin U\. 
For given values of ei, e 2 , fco (or u>) and and for an ar- 
bitrary function e(l), wc integrate the nonlinear ordinary 
differential equations (|10|l from / = to I = L numeri- 
cally, using the initial conditions The reflectivity TZ 
and the transmissivity T are obtained by 



TZ 



T 



e 2 - ei sin 



e2^/el cos 6* 



(12) 



In Models A and B, we do not need the equations for t 
and T, since the transmission is identically zero in those 
cases. In Model C, we note that the dielectric constants 
in regions I and II are identical. In that case, the initial 
conditions are simplified to r(0) = and t(0) = 1 and 
the transmissivity is given by T = \t\ 2 . 



If mode conversion occurs, the energy of the incident 
wave is absorbed into the inhomogeneous plasma, even 
when the damping parameter vanishes. In Models A and 
B, the mode conversion coefficient A, which measures the 
wave absorption, is obtained by A = 1 — TZ. In Model C, 
A is equal to 1 — TZ — T . 



The IIM can also be used in calculating the field ampli- 
tude B(z) inside the inhomogeneous plasma. We consider 
the B field as a function of both z and L: B = B(z, L). 
Then we obtain 



dB(zJ) U{1) 1 

— - — = ly/Tx &o cos 6 < - 

ol I ei 2 



e(0 
ei 



- 1 



1 777 tan 2 ( 



e(/) 



[l + r(l)}\B(z,l) 



(13) 



For a given z (0 < z < L), the field amplitude B(z, L) is 
obtained by integrating this equation from I = z to I = L 
using the initial condition B{z, z) = 1 + r(z). 



IV. RESULTS ON MODEL A: LINEAR 
DENSITY PROFILE IN ENTIRE SPACE 

The mode conversion in a warm, unmagnetized plasma 
with a linear density profile in entire space has been stud- 
ied analytically by Hinkel-Lipsker et al.~ If we restrict 
our attention to the cold plasma case with zero temper- 
ature, our Model A is identical to the model studied in 
Ref. 25. We introduce the dimcnsionlcss variables 




/, . .2/3 . 2 

(k n A) sin 



z = z/A, L = L/A, 1 = 1/ A, ( = k A, 



r(0) 



cos + i 



sin 



COS i 



■-n/T+ 



sin 



(14) 



and solve the invariant imbedding equation for r in 
Eq. I|1(J|) with the initial condition 



(15) 



where we have used the condition e 2 = — e\. 

In the limit where L — > oo and r\ — > 0, we find 
numerically that the mode conversion coefficient A (= 



1 - 



(fcoA) 2 ^ 3 sin 2 9. That this has to be the case can be 
easily seen by rewriting our wave equation using the 



is a universal function of the parameter Q a 



FIG. 1: Mode conversion coefficient A as a function of Q a = 
(/coA) 2 / 3 sin 2 6, for Model A. Our exact result is compared 
with the numerical result of Ref. 24 and the analytical result 
of Ref. 25. 



variables introduced in Ref. 25, z = fco 2 ^ 3 zj A 1 / 3 , L 



dB 



(z-L/2 + ifj-Q a ^ B = 0. 



d 2 B _ 1 

dz 2 z - L/2 + it] dz 

In Fig. 1, we show the mode conversion coefficient as a 
function of the variable Q a . In obtaining this result, we 
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(k A) 2,3 sin 2 e (z-z R )/A 



FIG. 2: (a) Mode conversion coefficient A for Model B as a 
function of the incident angle 9 for several values of £ = fcoA. 
(b) A as a function of Q a = (fc A) 2/;i sin 2 9 for Model B. As ( 
increases, the curves converge to the universal shape in Fig. 1. 



FIG. 3: Absolute value of the magnetic field, \B\, as a func- 
tion of (z — zr) /A for Model B, where zr (= L — A) is the 
coordinate of the mode conversion point. £ is equal to 12.5. 
(a) 9 = 5° (b)6 = sin' 1 0.4 w 23° (c) 9 = 45° 



have used a sufficiently large L (~ 100) and a sufficiently 
small i] (~ 10~ 8 ) and the convergence was perfect. We 
compare our result with the numerical result of Ref. 24 
and the analytical result of Ref. 25. We find that the 
agreement with Ref. 24 is very good, though we believe 
that the numerical accuracy of our result is better. 

V. RESULTS ON MODEL B: LINEAR DENSITY 
PROFILE IN HALF-SPACE 

In order to obtain the mode conversion coefficient in 
Model B, we solve the invariant imbedding equation for 
r in Eq. i|10[l with the initial condition 

(L ~ l) cos 8 + i\/L- cos 2 e 

r(0) = { (17) 

(L- I ) cos 8 -i\/L- cos 2 8 

where we have used the values t\ — 1 and 62 = 1 — 
L. In this case, the mode conversion coefficient is not 
a universal function of Q a even in the limit of L — > 00 
and rj —> 0, since the linear density profile persists only 
in half-space. In Fig. 2, we show our exact results for 
several values of £ = fcoA as functions of the incident 
angle 8 and the parameter Q a . As C increases, the curves 
converge to the universal shape in Fig. 1. 

We have also calculated the spatial dependence of the 
magnetic field amplitude B(z) using the invariant imbed- 



ding equation (|13|) . In Fig. 3, we plot the absolute value 
of the magnetic field amplitude, \B\, for ( — 12.5 and 
several values of the incident angle 8 as a function of 
(z — z^) /A where zr (= L — A) is the coordinate of the 
mode conversion point. After a suitable rescaling of the 
variables, our Fig. 3(b) agrees quite well with Fig. 2(f) 
in Ref. 24. 



VI. RESULTS ON MODEL C: PARABOLIC 
DENSITY PROFILE IN ENTIRE SPACE 

The mode conversion in a warm, unmagnetized plasma 
with a parabolic density profile in entire space has been 
studied analytically by Hinkel-Lipsker et al.w^^^ If we 
restrict our attention to the cold plasma case with zero 
temperature, our Model C is identical to the models stud- 
ied in Refs. 26—28. In order to obtain the mode conver- 
sion coefficient, we solve the invariant imbedding equa- 
tions for r and t, Eq. IjlOfl . with the initial conditions 
r(0) = and i(0) = 1. In the limit where L — > 00 
and r\ —> 0, we find numerically that the mode conver- 
sion coefficient A (= 1 — |r| 2 — \t\ 2 ) is a universal func- 
tion of the parameters Q c = fcoAsin 2 6>/(l + 5) 1 / 2 and 
A 2 = k A5/(l + 5) 1 / 2 . That this has to be the case can 
be seen by rewriting the wave equation using the vari- 
ables introduced in Ref. 28, z = k 1/2 (l + (5) 1 / 4 z/A 1 / 2 , 
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k A sin 2 e k Q A sin 2 e 



FIG. 4: Absorption coefficient A for Model C as a function of 
Q c = fcoAsin 2 9, for several values of the damping parameter 
r\ c . The wave frequency is fixed at the peak plasma frequency, 
therefore 8 is equal to zero. The value of L used in the calcu- 
lation is 400. Our exact result is compared with the analytical 
result of Ref. 27 represented by dashed lines. 
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FIG. 5: Reflectivity TZ for Model C as a function of Q c , for 
several values of r/ c . 8 is zero and L is 400. Our result is 
compared with that of Ref. 27 represented by dashed lines. 



L = k 1/2 (l + S^L/A 1 / 2 , Vc = k A V /(l + 5)V 2 : 
d 2 B 2(z- 1/2) dB 



dz 2 



(z-L/2) 2 -A 2 + i Vc 



dz 



+ 



(z-L/2) 2 -A 2 +i Vc -Q c B = 0. (18) 



When the wave frequency is equal to the local plasma 
frequency at the top of the parabolic density profile, cor- 
responding to S = 0, the quantity A vanishes unless 
i] c > and there is no true mode conversion. In this 
case, we will call A, which measures the absorption due 
to collisional damping, as the absorption coefficient in- 
stead of the mode conversion coefficient. 

In Figs. 4, 5 and 6, we show the absorption coefficient, 
the reflectivity and the transmissivity when S — as a 



FIG. 6: Transmissivity T for Model C as a function of Q c , 
for several values of r\ c . 8 is zero and L is 400. Our result is 
compared with that of Ref. 27 represented by dashed lines. 
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FIG. 7: Absorption coefficient when the angle of incidence is 
zero, A(9 — 0), for Model C as a function of the damping 
parameter r) c . 8 is zero and L is 400. 



function of the parameter Q c for several values of rj c . The 
scaled thickness of the plasma L used in the calculation is 
equal to 400. Our numerically exact results are compared 
with those in Ref. 27. The agreement is quite good when 
rj c = 0.0001, but becomes worse as r\ c increases. In fact, 
if we do the calculation for larger L values, we find our 
absorption coefficient gets bigger and the discrepancy be- 
tween our result and that of Ref. 27 increases gradually. 
This is reasonable because the damping occurs in the en- 
tire region of thickness L in our calculation, whereas sev- 
eral approximations limiting the effects of damping to a 
small region surrounding the peak point have been made 
in Ref. 27. We have calculated the absorption coefficient 
for the cases where the damping parameter r\ c is nonzero 
only in a narrow region of small width surrounding the 
point z = L/2. The discrepancy between the result of 
this calculation and that of Ref. 27 is much smaller than 
in Fig. 4, but there is no unique way of choosing a par- 
ticular width. 




FIG. 8: Absorption coefficient A for Model C as a function of 
the parameter A for several values of ij c . Q c is 0.05 and L is 
400. Our result is compared with that of Ref. 28 represented 
by dashed lines, (a) r) c = 10~ 8 (b) r) c = 0.001 (c) r) c = 0.01 



In Fig. 7, we plot the value of the absorption coefficient 
when the angle of incidence is zero as a function of the 
damping parameter. <5 is zero and L is 400. We find an 
almost linear increase of A(8 = 0), whereas it remains 
zero in Ref. 27. 

In Fig. 8, we show the absorption coefficient as a func- 
tion of the parameter A for several values of the damping 
parameter. L is 400 and Q c is 0.05. Our exact result 
is compared with the analytical result of Ref. 28. The 
agreement is quite good when -q c = 10~ 8 , but becomes 
worse as rj c increases for the same reason as in Figs. 4, 
5 and 6. If A is nonzero, the wave frequency is shifted 
from the peak plasma frequency and the mode conversion 

1/2 

occurs at two points, z — L/2 ± [5/ (5 + 1)] A, where 
S is related to Q c and A by S = (A 2 sin 2 0) /Q c . For a 
finite r) c , the quantity A includes the absorption due to 
both collisional damping and mode conversion. As i] c ap- 
proaches zero, A becomes the genuine mode conversion 
coefficient. The curve in Fig. 8(a) is a well-converged one 
which shows only the absorption due to mode conversion. 
The differences between Fig. 8(a) and Figs. 8(b) and 8(c) 
represent the absorption due to collisional damping. 

In Fig. 9, we plot the mode conversion coefficient, 
which is the absorption coefficient in the r\ c — > limit, 
as a function of A for several values of Q c . The values 
of L and r\ c used in the calculation are 400 and 10~ 8 
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FIG. 9: Mode conversion coefficient A for Model C as a func- 
tion of the parameter A for several values of Q c - tj c is 10 -8 
and L is 400. Our result is compared with that of Ref. 28 
represented by dashed lines, (a) Q c — 0.2 (b) Q c = 0.4 (c) 
Q c = 0.6 (d) Q c = 0.8 

respectively. We have verified that these results are well- 
converged in the sense that the results of the calculation 
for L values bigger than 400 and r\ c values smaller than 
10~ 8 arc indistinguishable from those shown in Fig. 9. 
We believe this is the first exact calculation of the mode 
conversion coefficient in the parabolic case. Our exact re- 
sult is compared with the analytical formulas of Ref. 28. 
Unlike the Q c — 0.05 case shown in Fig. 8(a), the dis- 
crepancy between our result and that of Ref. 28 is sizable. 
The shape of the curves is not very simple. For instance, 
we observe a double-peaked structure in Fig. 9(b). 

VII. CONCLUSION 

In this paper, we have presented a new version of the 
invariant imbedding theory of the wave propagation in 
stratified media and applied it to the mode conversion 
phenomena of p waves in cold, unmagnetized plasmas. 
We have obtained the mode conversion coefficient and the 
field distribution exactly for linear and parabolic plasma 
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density profiles for the first time. Using a recent general- 
ization of the invariant imbedding theory for the propa- 
gation of coupled waves in stratified media, we can apply 
our method to more general situations in a straightfor- 
ward manner.— In forthcoming papers, we will present 
the result of our study on the mode conversion in mag- 
netized plasmas and the temperature effects on both un- 
magnetized and magnetized plasmas. 
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